USCHEP/0311ib6 

CERN-CH/2003-285 

OHSTPY-HEP-T-03-014 



hep-th/0311264 



Improved Off- Shell Scattering Amplitudes 
in String Field Theory 
and New Computational Methods 

Itzhak Bars 

Department of Physics, USC, Los Angeles, CA 90089-0484 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

I.Y. Park 

Department of Physics, The Ohio State University, Columbus, OH 43210, USA 



Abstract 

We report on new results in Witten's cubic string field theory for the off-shell factor in the 4- 
tachyon amplitude that was not fully obtained explicitly before. This is achieved by completing the 
derivation of the Veneziano formula in the Moyal star formulation of Witten's string field theory 
(MSFT). We also demonstrate detailed agreement of MSFT with a number of on-shell and off-shell 
computations in other approaches to Witten's string field theory. We extend the techniques of 
computation in MSFT, and show that the j=0 representation of SL(2,R) generated by the Virasoro 
operators Lq,L±\ is a key structure in practical computations for generating numbers. We provide 
more insight into the Moyal structure that simplifies string field theory, and develop techniques 
that could be applied more generally, including nonperturbative processes. 
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I. INTRODUCTION 



In this paper we will report on new results in Witten's string field theory |l| for the off-shell 
4-tachyon amplitude that was not fully computed explicitly before. In previous computations 
the off-shell factor / (x) was obtained in an implicit form 2] ^ , while an explicit computation 
reported only the first two terms in an expansion in the vicinity of the integration limit 
of the parameter x 0. These results were obtained by using mainly conformal mapping 
techniques 0| and followed the methods of Giddings' original computation [^J of the on- 
shell Veneziano amplitude, which could be extended to off-shell under the guidance of the 
oscillator formulation of string field theory In this paper we will obtain a fairly 

comprehensive profile of the off-shell factor / (x) in its entire range < x < 1 by giving 
the explicit form to tenth order in the parameter x (first two terms in agreement with [^]), 
obtaining a plot in the full range, and computing the critical slope at a turning point in 
midrange that determines differentiability Contrary to our result, previously it was thought 
that the function / (x) was not differentiable at x = 1/2, where it was not well understood. 
These are achieved in the Moyal star formulation of Witten's string field theory (MSFT) 
U 



The usual approach of computation starts with a precise formulation of Witten's string 
field theory, such as the oscillator formulation, to derive a formal expression for a string 
Feynman diagram in terms of the cubic vertex defined in terms of the Neumann coefficients. 
After this step a jump is made to conformal maps from an analog model 0] and the real 
computation is performed by using conformal field theory, if the conformal map can be con- 
structed. The desired conformal maps can be found explicitly only in certain lucky cases, 
and the four point function is one of them. The conformal map procedure has been used 
virtually in every successful analytic computation, while the oscillator basis is directly pur- 
sued mainly with numerical studies using level truncation [18[ [19( because of the complexity 
of the Neumann coefficients. These have been some of the challenging features of Witten's 
string field theory in varions formulations 00 that, despite the beauty of Witten's bas.c 
action, have led to limited results in string field theory. 

Some of the complexities of other approaches are not present, or take an easier form 
in MSFT. In this formulation the string joining star product is the simple Moyal product, 
and this reproduces directly results in string field theory, in agreement with conformal field 
theory, but without mapping back to conformal field theory or other intermediate steps. 
Therefore some of the lingering problems in string field theory seem to be good testing 
grounds for MSFT. 

In this paper we apply the MSFT techniques to the off-shell 3-point and 4-point am- 
plitudes. We will derive the Veneziano amplitude directly from MSFT and determine the 
off-shell factor far more accurately compared to previous computations. In this process we 
demonstrate that MSFT agrees in detail with other approaches while bringing efficiency into 
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the computations, and also obtain new results in MSFT that other approaches could test. 
We view these computations as a preparation for our ultimate aim which is the investigation 
of nonperturbative string theory phenomena by using the simpler MSFT. We believe the 
techniques and insights developed in this paper will be useful for this purpose. In fact we 
find that some of the new results and insight gained by this work impact the way to compute 
nonpertubative quantities in string field theory as discussed in section VI. 

In the rest of this section we will describe briefly MSFT and introduce some notation. In 
section II we discuss off-shell 3-point functions. This is necessary to understand the behavior 
of the theory with respect to a cutoff in mode space, and to define the physical string coupling 
as opposed to the bare divergent coupling that appears in the action. The results we obtain 
for the off-shell 3-point amplitudes demonstrate detailed agreement between MSFT and 
conformal field theory or the oscillator formulation of string field theory. 

In section III we analyze the off-shell 4-tachyon amplitude. We show directly from MSFT, 
without connecting through conformal field theory, that on-shell we obtain the Veneziano 
amplitude. Furthermore we obtain the off-shell factor / (x) and compute the first ten terms 
in an expansion in powers of x, provide a plot in its full range < x < 1, and determine 
its differentiability at a turning point. Previously it was thought that the derivative was 
discontinuous at mid-range and only the lowest two terms were explicitly computed 0] (in 
agreement with our result). Our computation provides the most comprehensive information 
on the off-shell 4-tachyon amplitude produced so far in string field theory. 

In section IV we develop the details of the tools that allowed us to perform the com- 
putations in sections II and III and which would be applicable more generally to other 
computations. 

In section V we discuss the SL(2, R) generated by the Virasoro generators L , L±i, whose 
j (j + 1) = representation deeply underlies the structures that appear in our computations. 
By using some group theoretic properties of this very special representation of SL(2, R) we 
develop tools for computation in string field theory that are needed in our paper to generate 
numbers. We also compare the discrete and continuous Moyal bases which are simply two 
different bases of the j = representation. 

In section VI we discuss some of the impact that our present computations have on the 
non-perturbative landscape, and then conclude in section VII. In the Appendix we give 
further results produced through the techniques in section IV, and relations to Neumann 
coefficients. 

The action of Witten's cubic string field theory 0] in the MSFT formalism in the Siegel 
gauge is 

S(A) = - J d d x Tr C^A * (L„ - I) A + jA * A * Aj . (1.1) 
where Lq is given in Eqs. ()6.5l6.6|) . The zero mode ghosts have already been dealt with 
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|lq so they no longer appear. Hence the field A here is equivalent to the physical field 
in Witten's theory. The string field A(x,x e ,p e ) is written in a mixed position-momentum 
basis, which is equivalent to a Fourier transform of the purely position basis. In this basis 
string joining is represented by the usual Moyal star product [9] in a noncommutative phase 
space (x e ,p e ) labelled by the even modes e = 2,4,6, ■• • . The star product * is local in 
the midpoint variable x, and is independent for each e. This separation of variables and the 
simplicity of the Moyal star are the main conceptual and practical simplifications that lead 
to new progress by overcoming midpoint problems in other approaches and opening up easier 
computational techniques in MSFT. 

It has been shown in through some explicit computations that MSFT is in full 

agreement with other computational approaches to Witten's string field theory, including 
the oscillator formulation and conformal field theory 0-j2^. In particular, the Moyal star 
reproduces the Neumann coefficients that define the vertices in the oscillator formulation 
of string field theory 11 1 3]. Furthermore, the MSFT propagator has the usual free string 



spectrum. Therefore, even though it is a very different computational formalism, due to the 
one to one correspondence described in [3| we expect identical final results between MSFT 
and oscillator approach computations of any string Feynman graph. This expectation will 
be confirmed in detail in this paper. This demonstrates once again, and in greater detail, 
that MSFT is a precise representation of Witten's string field theory. 

The string in 26 dimensions is supplemented with two additional fermionic dimensions 
that describe the conformal ghosts b, c, with the appropriate generalization of the Moyal star 
product for fermions. The traditional perturbative string states (tachyon, vector, etc.) are 
identified through the usual expansion 

A(x cm ,x e ,p e ) = T(x cm )A (x e ,p e ) +V li (x cm ) (a£ 1 Ao(x ei p e )) + .... (1.2) 

where A (x e , p e ) is the perturbative vacuum string field configuration and at n is a differential 
operator representation of string oscillators in the space (x e ,p e ) A is a specific 

normalized Gaussian Tr(A -kA ) = 1, that represents the vacuum L A = 0, including 



The position basis is given in terms of the even and odd modes (x e ,x ) . The Fourier transform in the 
odd modes maps to the space Jx e ,Po) ■ Witten's string joining star product becomes a non-diagonal Moyal 
product in the space (x e ,p ) f9j. This is diagonalized by defining p e as an infinite combination of the p = 
p e T eo by introducing the special matrix T eo given in Ea.ljO)l. In this way we arrive to the noncommutative 
space (x e ,p e ) with a diagonal Moyal product whose meaning is string joining. We emphasize that this p e is 
a definition in terms of p 0l it should not be confused with the first quantized momentum that is canonical to 
x e , whose representation in this space is the derivative id Xe . However, interestingly this can be reproduced 
in the noncommuative geometr y (s tring joining) relation — id Xc A (x, x e ,p e ) = [p ei A (x, x e ,p e )]^ . A closely 
related continuous Moyal basis is obtained by orthogonality transformations, as will be discussed later 
in the paper. 
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ghosts [3] (see also fl^|). It is given by 

A (0 = |det4m | d/4 |det 4m x p 2/4 exp (-v^mo^C + ie^m^-^) , d = 26, 

(1.3) 

where £ m are the non-commutative coordinates written as a doublet for each e, with the 
bosonic part £f = (x^,p^) for matter and fermionic part = (x b e /\/2a', — y/2a'pl), £ 2 = 
(Xg/v^a 7 , \/2a'Pe) for the 6, c ghosts 2 . Each pair (x,p) satisfies standard commuta- 
tion/anticommutation rules under the star product. These can be written compactly as 
= r f v<J ij f° r matter and {^f 1 )^"} = —%E mn <Jij for ghosts, where the symbol 
a = —6<j 2 is the Pauli matrix oi in the doublet space multiplied with a noncommutativity 
parameter 9. We take 9 = 1 by a choice of units. The mo that appears in the vacuum state 



Aq (£) is a matrix in mode space determined by Lo^4o = as shown in 11 1. Although mo is 
a simple matrix, we will not need it explicitly in this paper. 

The perturbative particle fields T(x cm ), V^(x cm ), etc. are expressed as functions of the 
center of mass x cm . The star product * is local in the midpoint of the string x, not in 
the center of mass x cm . Therefore, before evaluating the interaction for any perturbative 
field T(x cm ), V^Xcm), etc. one must first write the center of mass in terms of the midpoint 
Xcm = x + w e x e} where w e = —\/2(—l) e ^ 2 . This is a crucial step in computations of the 
star products. The midpoint had been a source of numerous problems in the split string 
formalism, and the resolution first was given in the context of MSFT in ^l| with the simple 
prescription just described. 

In the quadratic term the star product plays no role, and could be removed, as is usual 
for the Moyal product. The coefficient of the quadratic term in the action is chosen such 
that the particle fields are correctly normalized (after taking into account the definitions of 
the trace Tr including ghosts and the Virasoro operator L as given in [l^ |) 

SauadraUc = ~ J d*X Q^T^T + + • • • ^ (1.4) 

= j d d x l - (d Td T - VT ■ VT) + ^d V»d V, + ■■■ (1.5) 

Then from the cubic term one finds the Feynman rules and computes the amplitudes as 
described in For each external line we insert a string field representative of a 

particle. In particular an incoming tachyon with momentum is represented by the string 
field A Q (x e ,p e ) e cm , where A is the normalized vacuum field given above, and e lk ' Xcm is 
the center of mass plane wave which is part of T (x cm ). To compute its interactions one must 
write it in the form A (x e ,p e ) e lk -( x + w ^ x ^) which is a shifted Gaussian in (x e ,p e ) space. These 



2 Relative to 0]Q we are improving notation by introducing £ m , with m = 1,2, and the Sp(2) metric 
£mn, for the sake of connecting with an upcoming paper in which we will discuss some useful hidden 
symmetries that connect matter and ghosts in the MSFT formalism, and further simplify the structure 
and computations. 
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details are fully explained in where it is also shown how to compute the star product 



and the trace with efficient methods based on a monoid algebra of shifted Gaussians 
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In the expressions below the constant matrix t eo and vectors w e , v Q in even/odd mode space 
(and related matrices T eo , R oe ) are fundamental matrices in MSFT that encode the joining 
of strings and are given explicitly in Eqs. (j4.3J4.4j) . The matrices k q , k s are diagonal 

matrices that represent the odd/even oscillator frequencies n Q = diag (o) , o — 1, 3, 5, • - • and 
K e = diag (e) , e = 2, 4, 6, ■ ■ • as in Eq. (j3.6J) below. A bar on top of a square or column matrix 
symbol, such as i, v, etc. means the transpose of the matrix. In certain computations, to 
avoid associativity anomalies these infinite matrices must be replaced with their regulated 
N x N version as given in 

0QQQ 

and then N must be sent to infinity at the end. The 



regulated matrices obey some nice algebraic properties which are also shared by the infinite 
matrices, thereby permitting analytic computation in the finite N version. A particular form 
of the regulator which we have found useful in some computations is given in footnote El 

It is necessary to use the regulator in those computations where we suspect anomalous 
behavior, but otherwise the unregulated matrices can be used, as we will do for most of the 
computations in this paper. It is interesting to note that the regulated matrices with only a 
few modes (small N) reproduce approximately most of the numerical results we obtain with 
more sophisticated methods at = oo. An example is the Neumann matrices as given in 
the appendix of 15], and many of the numbers computed in the current paper, although we 
do not make the effort to demonstrate this point in this paper. 



II. OFF-SHELL 3-POINT FUNCTIONS 



In this section first we briefly outline the off-shell 3-point functions to establish the relation 
between the bare coupling g that appears in the action of MSFT and the on-shell tachyon 
coupling g, which is identified with the string coupling. The relation between the two involves 
a factor which diverges with the number of modes 2N as go ~ (2N) 3 ^ 2 g as will be explained 
later. All other scattering amplitudes are proportional to some power of the bare coupling 
go, and this must be first written in terms of the finite on-shell tachyon coupling g. After 
this step, it is seen that all amplitudes are finite. In this process one finds a renormalization 
factor in front of the amplitudes multiplying a power of the on-shell coupling g. An example 
of this is the factor g 2 (27/16) 3 /4 for the 4-tachyon amplitude in Eq. (j3.5|) below. 

The off-shell 3-tachyon amplitude is obtained by using the Feynman rules (note x 3! = 
2go) and inserting the tachyon field for the external leg in the Feynman graph, leading to 
the expression 

9l23 (ki) = 2g J d d xTr (A e ikl - (s+WXe) * A e lk2 - (x+WXs) * A e ik3 < x+WXe) ) . (2.1) 
The x integral gives a momentum conservation delta function (27r) 26 5^ {k\ + k^ + k^), while 
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the trace and star products are easily evaluated by using the bose/fermi monoid rules devel- 
oped in 11| 3- ^ e resu ^ is 



/27\ hh a '( k l+ k i+ k i) 
9123 (h) = g ( —J , (2.2) 

multiplied with the delta function. This is in full agreement with previous studies of 
the off-shell 3-tachyon amplitude 00|[3]. It gives precisely g on shell by definition 
<?i23 {ki) | a 'fc?=i = 9- The remaining off-shell factor in our calculation initially has the 
form exp(o;(|-|a'(A;? + A| + A|))), with u = -L (w, 0) am (3 + mg)" 1 (w, 0) T . This 
is obtained easily by evaluatin g t he star products and trace in Eq. (j2.1|) by using the 
simple monoid methods in jll|jl5|. Using the definition of mo this simplifies to lu = 
2vko 1/2 (3 + ity 1 k 1 ^ 2 v, which is then evaluated as u = In (27/16) later in this paper in 
Eqs. (|4.13l4.2Uj) . to produce the result in Eq. (|2.2|) . Similarly, in the same computation, the 
relation between the on-shell g and the bare coupling g initially takes the following form 

g Q = |^(g)3/2 det( ^ )d/2 ^^3^^2/2 | detmo pd/4| detm -l|-2/4 ^ j t 

becomes 

go = I 9 ^V / \ detW) -<«»/Met(i±MV det fi±iMV 2 , d = 2e , ( 2 .3) 



2 a \16J v v " V 4 / V 4 

The right hand side of Eq. (j2.3|) is divergent as (2iV) ( ' d ~ 6 ' )//8_d ^ 18+4 ^ 9 which becomes (2A r ) 3 ^ 2 for 
d = 26 as will be later shown in Eq. (|4.32|) . However this divergence is everywhere reabsorbed 
into the definition of the on shell g just as in the 3-tachyon case. After this step, there still 
remains similar determinants that individually produce a divergence or zeroes at large N, 
but they combine together to give finite answers magically as long as d = 26. We will see an 
example of this impressive fact below in the 4-tachyon amplitude (see the determinants in 
Eq. ()3.5|) and the computation in Eqs. (|4.241 )). 

A similar computation can be performed for the vector particle for an incoming wave with 
momentum which is represented by the string field 

A v (x cm , x e , p e , k) = e ik - x ™A (x e , p e ) (k) (p»T) 1 v 7 ^ 7 , (2.4) 

where e fl (k)is the polarization and the last factor p^T el \/2a' resulted from applying the 
oscillator U| 1^ a / i 1 A (x e ,p e ). For example, the coupling between the vector and two 
tachyons is (omitting the Chan-Paton factors at the string ends) 

g 12Vs (h) = 2g J d d xTr (^A e ik ^ x+wx ^ * A e ik2 < x+WXe) ★ A e ik ^ x+wx ^ (k 3 ) (j^T\ v 7 ^ 7 ) . 

(2-5) 

Following the monoid methods in ll| l^J that led to Eq. (j2.2j) this is evaluated in almost the 



same way, giving 



1 /27\ 1 -h a '( k i+ k 2+ k l) 



gi2v 3 (h) = -g(h - k 2 ) -e(k 3 ) f — ) (2.6) 
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with an implied momentum conservation delta function as in Eq. (J2.2j) . This result is in full 
agreement with the on-shell result in (for a'kf = a'k^ = 1 and a'k^ = 0), as well as 
the off-shell results summarized in [2]. To arrive at this result we needed to perform the 
computation of the following numerical coefficient 

^^^^-((^M-K- (27) 

This computation, and those for similar quantities will be given later in this paper in 

We see that MSFT is in full agreement with previous computations in string field theory 
for off-shell 3-point functions. Previous methods used conformal maps .2] and conformal field 
theory, while MSFT uses the Moyal product, and produces the same results with considerably 
simpler methods. One can go on computing very simply off-shell 3-point couplings for any 
other perturbative or non-perturbative fields by similar MSFT techniques. 



III. OFF-SHELL 4-TACHYON SCATTERING 



One of the new results in our paper is a proof that the MSFT 4-tachyon amplitude 
produced by the Moyal product does give directly the Veneziano amplitude, without using 
conformal mapping techniques. A second new result is a comprehensive understanding of 
the off-shell factor by going considerably further in the parametric expansion, obtaining a 
plot of the function in the full range, and determining its differentiability at a turning point. 

In MSFT we have been developing analytic methods of computation in string field theory 
by using directly the Moyal star product. Both the oscillator formalism and MSFT provide 
an expression for any off-shell amplitude, including loops. Although the starting point and 
intermediate steps a re quite different, it has been argued that generally we expect agreement 
in the result 



step s ar< 

00. 



The computation of the 4-tachyon off-shell amplitude in MSFT was performed in [lc 
There are 8 diagrams that correspond to various permutations of the four external legs. The 
s-channel diagram \ > — <±is denoted as 12^34- The mathematical expression for this ampli- 
tude was obtained in 0|0|. The t-channel diagram is given by the cyclic permutation of 
the external legs 41A23, and its mathematical expression amounts to exchanging s,t in the 
previous result. The remaining permutations that do not change the s,t channel properties 
are denoted as 01 A 43, 43A21, 34^12 and 14A32, 23 An, 32 A4 respectively. After a brief com- 
putation (as in 13]) it can be shown that these give the same amplitudes as the initial s,t 
diagrams. Therefore the sum of the diagrams for the s, t channels produces a factor of 4 

2 V 3\ „oo 

A (s, t) = 4 ( \> - <\ + 1 ]=4/ dre T ( 12 A U (r) + 41A23 (r)) 

1A4 
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CO 



dr e - 2a W +4 ^Wa 4 (r) e (*)+W)EtiK-i) (3.1) 

x e ( a 's+2)(r+a(r)+2/3(r)) e (a't+2)2/3(r) + / ^ j\" _ 

where s,t are the Mandelstam variables 3 s = — (ki + k^) 1 and f = — {k\ + kj )\ The func- 
tions a (r) , /? (r) , 7 (r) , a 4 (r) initially given in |l3| were further simplified in 3]. Here, as 
a first step, we simplify these functions further with a few algebraic steps and put them into 
the following form which will be convenient in our analysis 

a (r) = u - be~ TK (l - Me"™) ^ b, (3.2) 
(3 ( r ) = bP {o) e- TK (1 - Me-^Me'^y 1 P (o) 6, (3.3) 

7 (r) = — w - bP {o) e- TK (1 - Me' TK ) " X P (o) 6, (3.4) 



-TK 



1 /27\ 3 det f 1 ~ Me~ TK Me~ 

a 4 ( T )=A-g 2 — ^ ^r, d = 26. (3.5) 

41 ; 4 y Vie; [det (1 - Me— Me— )] d/2 1 ' 

The overall factor of 4 in 04 is the factor in Eq. (J3.1|) . In string mode space labelled by even 
and odd positive integers e = 2, 4, 6, • • ■ and o = 1, 3, 5, ■ ■ • , we have defined the matrices 
k, M, M, p(°\ vector b and scalar cu as follows (a bar above a symbol means matrix transpose) 

/. = 1 1. m= I v:i -" ; "' ., 1. m= j 1 (3 . 6) 





\ 1+3U 00' 



(3.7) 



Here the diagonal matrix k with integer eigenvalues n n = n, represents the spectrum of string 
oscillation frequencies (eigenvalues of Lq), P^ is a projector onto the odd modes only. The 
other quantities are all constructed from the matrix t eo and vector v Q which themselves are 
built from the frequencies n n as given in fo r any N. In the large N limit they take 

the form in Eq. (j4.3|) . The combination of t eo in the form of the matrix M arises frequently 
in the interactions in the matter sector, while the matrix M occurs in the ghost sector 4 . The 
even and odd vectors b e , b Q are given by 

2 ^o 1/2 A ,b e =( t-^-^v) , (3.8) 



3 + tt J ' V 3 + tt 



3 The Mandelstam t should not ne confused with the matrix t eo . 

4 The matrices M, M is are simplified forms of -M {0) , -CX^ where M {0) , were identified in 

as some of the Neumann matrices for the 3-point vertex in the matter and ghost sectors respectively. 
Furthermore, it was shown that all Neumann matrices, for all n-point vertices, are explicit functions of 
the matrix t eo as obtained in jll| |l "| . Hence t is the fundamental matrix that determines all interactions 
in string theory. 
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where k is the odd part of the matrix k. Finally u, which is the t —>■ oo limit of a (r) is 
simplified to the form 

u = b (1 - M)- 1 b = vk^j^k^v (3.9) 
by using the definitions of b, M given above. 

From these quantities we compute the functions a (r) , (3 (r) , 7 (r) , a 4 (r) which in turn 
determine the 4-tachyon scattering amplitude off-shell as well as on shell. We will show that 
the on-shell amplitude for a'kf = 1 reduces to the Veneziano amplitude given by the beta 
function 

, , 2 T(-a's-l)T(-a't- 1) 
A shell s, t) = g 2 1 , J \ 3.10 
T (— a's — a't — 2) 

On the other hand the off-shell expression above goes beyond conformal field theory which 
only gives information for on-shell strings. 

The aim in the rest of this section is to show that the off-shell 4-tachyon amplitude in 
Eq. (|3.1|) can be rewritten in the following form and then compute the function / (x) 



00 



A(s,t) = -g 2 ^ dT^-(f(x(r)))^^- 1 ) (3.11) 

(x(r)r' s - 2 (l-x(r)r' 4 - 2 + ( S -t) 

= g 2 f dx x~ a ' s - 2 (1 - x)- alt - 2 (f (x))^^?- 1 ) . (3.12) 
Jo 

After the change of integration variables from r to x in the form of Eq. (j3.12j) we see that 
the off-shell amplitude is consistent with the on-shell Veneziano amplitude (beta function in 
Eq- flTIDD ) when a'kf = 1. 

We will show that the change of variables from r to x is such that x (0) = 1/2 and 

x (00) = 0. Then we see that the s-channel amplitude 12^.34 (t) —\> <4 contributes to 

the range < x < 1/2 while the t-channel amplitude 41A23 (r) contributes to the range 
1/2 < x < 1 after a change of variables x — > (1 — x) . For the first form in Eq. (|3.1j) to 
agree with the second form in Eq. (j3.11|) it is required that a (r) , (3 (r) , 7 (r) , 04 (r) conspire 
to have remarkable relations among themselves so that they can be written as functions of 
the same x (r) . Thus we need to prove that the following relations are satisfied (which also 
define x (r) in terms of a (r) , (3 (r) , 7 (r) , 04 (r)) 

e -T e -a(T)-2/»(T) = j _ e -2/3(r) = £ ( r ) ? (3.13) 
rfX ( T ) = ° 4 ( T ) c -r-2a(r)+4 7 (r) /g 

Note that 04 (r) depends on ti = 26 while the other quantities are independent of the number 
of dimensions. Hence if the relation holds for d — 26 it cannot hold for other dimensions. 
We will prove below that the relations in Eqs. ()3.13I3.14|) are indeed true. This makes d = 26 
unique. 
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Once the relations are proven, then we learn that the off-shell factor / (x) is given by 

f(x(r)) = e 7{T)+2/3{r) . (3.15) 

To write / (x) in terms of only the parameter x in the integral representation of the off-shell 
amplitude in Eq. (j3.12|) . the relation between x,r given in Eq. ()3.13|) needs to be inverted 
r = t(x). We will perform the inversion and will construct the function / (x) as a series 
expansion in powers of x. It turns out that a few terms in the expansion already give the 
necessary information to obtain a sufficiently accurate representation of the function / (x) , 
and hence of the full off-shell 4-tachyon amplitude. 

Let us first prove that the relations in Eqs. (j3.13j) hold at the integration limits r = 0, oo. In 
the next section we show how to compute the functions a (t) , /3 (r) , 7 (r) , a 4 (r) at r = 0, 00. 
We find in particular 

ot(0) = 0, (3(0) = -ln2, 7(0) = — ln2, (3.16) 

27 ! 16 ( 27\ ^ 

a(oo) = ln— , /5(oo) = 0, 7(00) = -In— , a 4 (oo) =^ 2 ( — j . (3.17) 

From this we see that indeed Eqs. ()3.13|) is satisfied at both limits r = 0, 00, and we also 
determine 

x (0) = -, x(oo) = 0. (3.18) 

This shows that the s-channel amplitude 12^34 [j) is associated with the range < x < 1/2 
and the t-channel amplitude 41A23 (r) contributes to the range 1/2 < x < 1 after the change 
of variable x — > (1 — x) . 

Next we examine the relations for more general values of r. The form of the equations 
in Eqs. (j3.2H3.5J) was developed to examine an expansion in powers of e~ T . In the next 
section we show how to compute the coefficients for the large r expansion of the functions 
a (r) , (3 (r) , 7 (r) , a 4 (r) . We find the following analytic result 



. 16 8 _ T 2 2 19 _ 2t 2 5 7 2 _ 

a (t) = - In — -e 6 ' 

K J 27 27 (27) 2 3(27) 3 

/ n 1 1 16 8 _ T 40 _ 2t 2 5 7 2 _ 3r 

7 ( r ) = - In e ? e d 

n 1 2 27 27 (27) 2 3 (27) 3 

2 3 829 ^ 2 4 167 x 229 



-e 



-4r 



(27) 4 5 (27) 



5 



e~ 5r + O (e^ T ) , (3.21) 



M^^g) 3 ( 1 + ^e- + ™e- + OM ). ,3,2) 
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To obtain the expansion to this order it is sufficient to compute the coefficients 
bi, &2, b 3 , 6 4 , 65, M11, M13, M 2 2, M11, M3, M 2 2 defined above. These results were obtained an- 
alytically without much effort. Our analytic results above to order O (e~ 6T ) are already quite 
adequate to construct / (x) . It is possible to easily extend the expansion by inserting the 
results for the higher coefficients provided in the appendix into Eqs. (j3.2H3.5j) . In fact we 
have constructed an algebraic computer program that does this, and using it we have double 
checked our analytic results above and extended our computation to higher orders. We will 
report on some of the higher order results below. 

The relations (|3.13l3.14j) can now be verified directly by inserting the large r expansions 
for a (t) , ft (r) , 7 (r) , a 4 (r) given above, and re-expanding in powers of e~ r up to order 
O (e~ 7r ). From either the first or second term in Eq. ()3.13|) we obtain the same expression 
for x (r) , namely 

_ 128 -2t _|_ 64 -3t _ 19 456 -At 
T (-r\ = I 27 729 ' 729 531 441 e | (o oq\ 

X V ) I , 387296 -5t 1733120 -6r , f) ( r ~7r\ l> I •>•-■>; 

14 348 907 129140163 \ 



59. 259e" T - 17. 558e" 2r + 8. 779 le~ 3r - 3. 661 0e" 4T 
+2. 699 le" 5T - 1. 342e~ 6r + O (e~ 7r ) 



x 10" 2 (3.24) 



Furthermore, the right hand side of Eq. (J3.14|) gives precisely the derivative d T x (r) of the 
expansion in Eq. (J3.23|) 

n ( T \ ( 16 -r I 256 -2t 64 -3t i 77824 -4r\ 

" 4 V / -t-2«(t)+47(t) 27 e ^ 729 e 243 e ^ 531 441 e | ocN 

n 2 e I 1936 480 -5r i 3466240 -6r i n ( r~ 7r \ /' 

i>T \ 14348907 C " r 43046721 e V > J 

This proves that a (r) , ft (r) satisfy the relation in Eq. (j3.13|) . and a (r) , 7 (r) , a 4 (r) satisfy 
the relation in Eq. (j3.14|) . at least up to order O (e~ 7r ) . Actually, as already mentioned, with 
a computer program we have shown that the relations hold to much higher orders. These 
results are convincing that a (r) , ft (r) , 7 (r) , a 4 (r) are all expressed in terms of the same 
function x (r) . We emphasize that the number of dimensions d which appears in a 4 (r) in 
Eq. (|3.5|) must be d = 26 to satisfy the relations. 

Note that at r = we expect x (0) = 1/2 exactly as in Eq. (j3.18|) . and indeed by replacing 
e~ r — > 1 the expansion (|3.23|) we obtain x (0) = 0.481 76 which implies that the expansion 
captures an accurate representation of the full function. Note from the trend in Eq. ()3.24|) that 
the next O (e _7r ) correction will bring the value much closer to the exact answer x (0) = 1/2. 
With the computer program we have shown that this value becomes x (0) = 0.492 95 by 
computing x (r) to order O (e _10r ). 

Finally we compute the off-shell factor, which is given by 

/ (x (r)) = e^r) = ^ (1 " a:(r)) (3.26) 

The expression on the right hand side defines the function C which is given here for com- 
parison to the old literature. Thus, by substituting 1 — x (r) = exp (—2/5 (r)) we can write 
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C (1 — x (t)) = 2e 7 ( T ) + ^ T ). First we obtain the expansion of C (1 — x (r)) by inserting our 
computation of 7 (r) and /3 (r) . The result is 



C (1 - , (r» = > = g) ^ (. - 1 (») 2 e - - 1 (g 2 + O (e-) 

(3.27) 

Next, by using Eq. ()3.23|) we rewrite this result in terms of x (r) so that the function 
C (1 — x (t)) is given in terms of x. The result is 



C(l-x(r))-( A X /2 ( l-|(^(-)) 2 -|(^(-)) 3 

1 [ )) W I -H(^(r)) 4 -|i(x(r)) 5 + 0((x(r)) 6 ) 



(3.28) 



1. 539 6 - 0.240 56x 2 - 0.240 56x 3 
-0.234 74x 4 - 0.228 91x 5 + O (x 6N 

These coefficients are determined by demanding that the expansion in powers of e~ T matches 
the one of 2e 7(r)+/3(r) . In this form our result for C is in agreement with what is found in the 



old literature |3[ where C was computed up to the second term — J^x 2 by using Mandelstam's 
conformal mapping techniques. We, of course, used the very different Moyal star technique 
and obtained the same result, but we also easily went further by obtaining the higher order 
terms, which is a new result given in this paper. As mentioned, with our technique it is very 
easy to compute to even higher orders (see Eq. (j3.33J) ). 

Since the range for xis0<x<l/2 the expansion given above for C has good convergence, 
so we expect that we have obtained an accurate representation of the full function C in the 
relevant range. As a test let us compare the exact value of C at r = (or at x — 1/2) which 
we can compute exactly C(l - x (0)) = 2e 7(0)+/3(0) = V? = 1. 414 2 after using Ea. (j3~ToT) . By 
evaluating the expansion above at x — 1/2, we obtain C — 1. 4276, which confirms that the 
expansion does capture the function almost fully in the entire range. 

Let us now turn to the full off-shell factor / (x) . The exact expansion of this function in 
powers of x up to O (x 6 ) becomes 

n Cl — t\ A /l xi T 4- J-t 2 -A- _ i29_^,4 \ 

y> 2(l-xf 2 3v/3l -^x 5 + 0(x 6 ) ) y - ' 



1 + 0.5x + 0.218 75x 2 + 0.07813x 
3^ I -0.01575x 4 - 0.08624x 5 + O (x 6 ] 



(3.30) 



At the end of the range x — 1/2 (i.e. r = 0) the exact value of this function should be 1, 
since 

/ (1/2) = f(x (0)) = e 7(0)+2/3(0) = 1, (3.31) 

where we used Eq. (j3.16|) . The series approximation in powers of x given above to order 
O (x 5 ) produces / (1/2) = 1. 0090, which is better than 1% accuracy. 
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To get a better feel of the function / (x) we plot the expansion in Eq. (J3.29|) for the 
appropriate range for x which is < x < 1/2. This is shown by the solid line in Fig.l. 

l: 
0.95: 
0.9: 
0.85: 

0.8: 

0.1 0.2 x 0.3 0.4 0.5 

Fig.l 

The figure suggests that, for the relevant range, as a guide for the eye we can compare / (x) 
roughly to a linear function / (x) 

where the slope is chosen to guarantee the exact values at both ends of the range, namely 
at x = and 1 at x = 1/2. This case corresponds to the dotted line in the figure. The 
plots suggest that the exact curve / (x) lies somewhere close to the solid line and the dotted 
line, and that in any case we have obtained a fairly good approximation in the entire range 
with the series expansion in Eq. (jri.29j) . 

By combining the s and t channel contributions we extend the range of integration to 
< x < 1 as explained above. Thus, in the range 1/2 < x < 1 the function / (x) is given by 
substituting a; by (1 — x) in the expression of / (x) given above. Namely, the plot of / (x) 
in the full range < x < 1 is given in Fig. 2. 

i: 

0.95 
0.9 
0.85 
0.8 

0.2 0.4 x 0.6 0.8 1 

Fig.3 

The exact function is expected to lie somewhere close to the solid and the dotted lines. 
The shape of the curve begs the question of whether the exact / (x) ever crosses the dotted 
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line. To provide an answer we used our algebraic computer program in which we plugged 
in the contents of the Appendix into Eqs. (j3.2H3.5|) and obtained the following higher order 
expansion (we do not give the details of the expansion for a(r) , (5 (r) , 7 (r) , a 4 (r)) 

4 / 1 4- It 4- -Lt 2 4- A r 3 _ J^9_ T 4: _ 1413 5 _ 40 973 6 \ 

f ( T \ — [ 2 32^ ~ 64 8192^ 16384" t ' 262144^ | (o oo\ 

o /q I 124459 7 186 841 777 8 547864633 9 , q / 1Q\ I" ^o.oo; 

°V<J \ 524288 536870912^ 1073741824 x ~ r ^ v x / / 

The plot of this more accurate expansion lies slightly below the solid line in the figure, and is 
extremely close to it in the regions that are not near to x — 1/2. It does cross the dotted line, 
but it does it at mirror points closer to x = 1/2 compared to Fig. 2, and finally reaches the 
value / (1/2) = 1. 003 9 at x = 1/2. So the new corrected peak is between the previous value 
of 1.0090 and the exact value 1.0090 > 1.003 9 > 1. This analysis is consistent with the 
possibility that the exact function / (x) may cross the dotted line somewhere close to x = 1/2 
before settling into the value / (1/2) = 1. This issue of crossing or not crossing is relevant for 
determining whether the exact function is different iable at x — 1/2. If it never crosses the 
dotted line then it must make a cusp and be non-differentiable. Such a non-differentiable 
function may be a peculiarity of Witten's theory in which a particular conformal gauge has 
been effectively fixed by distinguishing the midpoint. We remind the reader that off-shell 
amplitudes generally are not gauge invariant. 

It seems therefore interesting to investigate the slope of / (x) near x — 1/2. A glimpse of 
the slope on the left side of the possible cusp is obtained by computing the derivative of the 
expansion above and evaluating it at x = 1/2. We obtain d x f (x) \x=i/2 — 0.49849 consistent 
with a cusp. However it is not clear how much we can trust this number, because if we 
examine the size of the coefficients in the expansion of the derivative 

/ 38. 49 + 16. 839 (2x) + 4. 510 5 (2x f - 0.606 11 (2x) 3 - 2. 074 7 (2x) 4 \ 1Q _ 2 
^ ~ \ v -2.256 (2x) 5 - 1.998 7(2x) 6 - 1.6744(2x) 7 - 1.380 9 (2x) 8 + O(x 9 ) J X 

(3.34) 

we see that convergence is not very fast, and the exact value may turn out to be quite 
different at 2x = 1 since there are an infinite number of terms to be summed. Note that the 
tendency of the neglected terms is to be negative thereby reducing the value of the slope 
given above. To settle the question we must examine an expansion near r = rather than 
near r = oo, but the reason we did not do this so far in this paper is because the expansion 
does not seem to exist due to divergent derivatives near r = 0. However, precisely these 
divergences turn out to settle the issue as follows. 

We can compute the derivative (d x f {x)) x=1 / 2 from d T f (x (r)) | T= o as follows 

* (0) (dj (x)) x=1/2 = d T f (x (r)) | r=0 = ( 7 (0) + 2$ (0)) e^+^°\ (3.35) 

where we have used the chain rule on the left hand side and Eq. (j3.26)l on the right hand 
side. From Eq. (|3.13J) we have x (0) = 2(3 (0) e~ 2/3 (°\ Inserting the exact values for j (0) , j3 (0) 



16 



given in Eq. (j3.16j) these equations become 



x(0)=$ (0) , x (0) (d x f (x)) x=1/2 = 7 (0) + 2(5 (0) . (3.36) 

By differentiating Eqs. ()3.3l3.4|) we obtain expressions for 7 (0) , (0) 

/3(0) = -b (1 - M fl 2 ) -1 (k c + M k M ) (1 - Mq) 1 b a , (3.37) 
7 (0) = b Q (1 - M^- 1 Ko (1 - Mo)' 1 b , (3.38) 

where M Q = |^ is the odd part of the matrix M. After srtaightforward algebra, we find 



and x (0) , /3 (0) related to the same expression by 

i(0) = /5(0) = -i 7 (0)-i (3.40) 

We used vv = 1 which multiplied the last constant term —1/2. In passing we mention 
that we have also computed a (0) and found a (0) = 27 (0) after a little algebra. Through 
Eqs. (l4.16l4.17|) and Eqs. (j5.36H5.4U|) in the next section we show that the expression above for 
7 (0) is divergent at large N. So 7 (0) + 2(5 (0) = —1 is finite, but x (0) diverges. Therefore 
from Eq. (j3.36|) we find that the slope of / (x) at x — 1/2 is exactly zero 

(d x f(x)) x=1/2 = 0. (3.41) 

So, after all the off-shell function f (x) is continuous and differentiable at x = 1/2. With 
this result on the slope, the expansion in Eq. (j3.33|) . and the plot in Fig. 2, we have basically 
understood the function / (x) . 

This discussion provides the most comprehensive result for the off-shell scattering ampli- 
tude of four tachyons produced so far in string field theory. 



IV. COMPUTATIONAL TECHNIQUES 

In this paper we need to evaluate the quantities u, b D , b e , c e , c Q , M ee r, M 0o t, M ee >, M 00 i, and 
various determinants, as defined in the previous sections. These are examples of more 
general computations that come up in MSFT which are generically of the type 

(F (it)) 00 > , (F (tt))ee> , (tF (it)) eo , («) q , (tF (it) K^v)^ (4.1) 

det (F (it)) , vk~ 1/2 F (it) n~ l/2 v, vk~ 1/2 F (it) k F (it) k^v, etc. (4.2) 

where F (z) can be any function of the matrix it or it. In the present paper the functions 
F of interest are M = and M = and other simple ones constructed from them, 
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such as = 1 — M, (1 — M 2 ) , etc. As seen in the full r-dependent a(r) ,{3 (r) etc. 
there are more involved combinations of F (it) and the matrix e _TK which can appear in 
explicit calculations. The techniques discussed below do not directly apply to compute such 
quantities fully, but do apply to compute the expansion coefficients analytically and exactly 
when these quantities are expanded in powers of e" T . This is how we made progress in this 
paper. 

In our MSFT approach such quantities are initially defined in terms of the regulated 
vector v a and the regulated matrix t eo . By inserting the values of v a , t eo for any number of 
modes 2N, all of these quantities take on explicit numerical values. Generally it is seen that 
a few modes already give an answer pretty close to the N — > oo limit, therefore for a quick 
estimate of non-divergent quantities it is generally sufficient to use finite matrices for just a 
few modes N. With this approach we can evaluate any quantity approximately fairly easily. 

The large N limit is sometimes subtle because of anomalies. Generally anomalies arise 
in the form of a zero multiplied by an infinity that comes from an infinite sum over matrix 
indices (the modes) producing x oo, so that there is a subtle finite contribution. If one 
is not careful the zero is first evaluated and the infinite sum that is later performed still 
gives zero, thus missing the finite contribution. To avoid anomaly subtleties in the N — ► oo 
limit we use the regulated version of v ,t eo in all computations (see [3]), and take the 
N — > oo only at the end of the computation. The regulated matrices v Q , t eo satisfy certain 
algebraic relations, also shared by the unregulated matrices, that make it possible to perform 
certain computations analytically. It has been demonstrated that this is a correct regulator 
consistent with other computational techniques in string theory. 

In certain computations, or parts of it, there are no anomalies, and it becomes possible 
to work with the infinite unregulated matrices v e , t eo directly at N = oo. Then it is possible 
to make a transformation to the basis that diagonalizes the infinite unregulated matrix t eo 
and do the computations in the N = oo diagonal basis. At N = oo the eigenvalues of 
t eo are continuous and are parameterized by a continuous variable k, and furthermore the 
transformation is known explicitly As we will discuss in detail in the next section, 

at iV = oo the discrete and continuous bases are just different bases of SL(2, R) in the 
j (j + 1) = unitary representation, in which t eo can be understood as the matrix elements 
of the operator t = tanh (L\ + in the discrete basis that diagonalizes Lq. The 

continuous basis is the basis in which this operator is diagonal. 

A priori it is hard to know when to expect an anomaly if one works directly at iV = oo. 
Therefore the regulator provided in MSFT is indispensable to insure that such subtleties 
will not spoil a computation. So one must be wary when using the continuous basis directly 
at iV = oo. We emphasize that the discrete Moyal basis with a regulator is the safe way to 
proceed in general, but one can make a transformation to the continuous basis for specific 
computations to evaluate infinite sums in parts of computations. This is the sense in which 
we use the continuous basis in the next subsection, and indeed we will find good use for it 



18 



in some of the following computations. 



A. Computations via transformation to continuous basis 

In the following computations we are directly at iV = oo and work with the unregu- 
lated fundamental matrices T, R, v, w in MSFT that encode string joining (for a particular 
regularization see footnote [HJ) 

w e = V2 (t)- e+2 , v = — ^ , t eo = J~eT eo ^ (4.3) 

7T O Jo 



7T(e 2 — o 2 )' 7ro(e 2 — o 2 ^ 



As explained in ll| 12j, the matrix t is diagonalized £ eo = (V e rVoj eo by the orthogonal 
matrices V Q ,V e that act on the odd/even sides (a bar means transpose). The eigenspace is 
labelled by k with eigenvalues 5 t k . The label k is continuous (can happen only for infinite 
matrices); it is in the range < k < oo, and the eigenvalues are given by ^ t k = 
tanh(7r/€/4) . The orthogonal matrices have matrix elements which are functions of k and 
the index o, e, namely (V ) 0K = V Q («) , and (V e ) eK = V e («) . Thus, the matrix equation 
teo = (V e ) eK r K (y o ) KO is written as 

/■oo 

t eo = / d«y e («) tanh (7T«/4) K («) (4.5) 



o 



The functions V e (k) , V a (k) satisfy orthogonality relations that correspond to the orthogonal- 
ity conditions on the matrices V e , V D . These functions are given explicitly by the generating 
functions in Eq. (|5.Hjl . In the next section we will clarify the role of these functions as the 
overlaps < e\n >= V e (k) /a/2 and < o\k >= V a (k) /a/2 between the states of two different 
bases of the j = representation of SL(2, R) generated by the Virasoro generators L , L ±1 . 
The first basis is the familiar one |n > labelled by the even or odd eigenvalues n = (e, 6) of 
Lq. The second basis \k > diagonalizes (L + i + L_x) \k >= k\k > . 

By evaluating the generating functions in Eq. (|5.3J) at k = and expanding in powers of 
z, one can see that V e (0) = 0, while V Q (0) is finite and related to v as V a (0) = \ Jlwv . 
Thus v is directly related to the state \k = >. We have evaluated the non-trivial infinite 
sum J2 ^/f^ m ^ ne nex ^ section. We give it below along with a few of the V e (k) , V (k) 
that are needed in our computations 

Vl f K ) = V 2 (k) = ; K " , V 3 («) = ~ /t + ^ 3 , (4.6) 

1V ; ysinhf 2y ' ^Ksinhf V; ^ sinh f V ; 



5 This continuous parameter k should not be confused with the diagonal matrix k of Eq. l|3.6|l . As much 
as possible we are trying to keep the notations that were introduced in different papers by independent 
authors. We hope the reader will discern them in the proper contexts. 
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Vl(k) 



4k sinh ^ttk 



V r (k) 



k (1 - §k 2 + -^k 4 ) ^K(k 



5k sinh ^ttk 



tanh 22£ 

(4.7) 

«; cosh 2 ^ v ; 



With this preparation we are ready to compute the desired generic quantities. Applying 
the diagonalization of t eo we have just described we can write t = V e rV , it = V t 2 V and 
ti = V e r 2 V e and express the desired quantities as integrals 



(F(it)) oa , 
(F(tt)) ee , 

(tF(tt))eo 

{F(it)K- 1 / 2 v) o 
(tF (it) «-V2 v ) 



dnV (k) F ( tanh 2 (^)) («) 

f//, r: (//) r ( tanh 2 v e , («) 

JkI; f/v)tauli ('™) F (tanh 2 V (k) 



GkV (k) F (tanh 2 



tanh^f 
k cosh 2 ™ 



f//,r: (/■,') tanh [ ^ F (tanh 2 



tanh^ 
k cosh 2 



UK 



/ 2 , 



(iKi 



tanh^ 



Kcosh 2 ^ 



1 Fftanh 2 ^— V 



det F (tt) = exp / cfKp (k) In 



V 4 A 



tanh^ 
Kcosh 2 ^' 



In the infinite N limit we have det (F (ti)) = det (F (tt)) . The quantity 

ZK 



p{K) 



(4.8) 
(4.9) 
(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 

(4.15) 



which appears in the evaluation of determinants will be derived in the next section in 
Eq. (j5.58|) . Here if> (z) = d z lnT (z) is the logarithmic derivative of the gamma function, 
and A depends on the regulator. It could be chosen as zero if we compare to a particular 
regulator in the discrete basis of MSFT. The same result has been obtained in a very differ- 
ent type of calculation in 00- The term that contains the number of modes 2N is the 
leading term independent of k. In physical computations, as N goes to infinity this factor 
cancels magically among matter and ghost determinants as long as d — 26, and the bare 
coupling go is rewritten in terms of the on-shell tachyon coupling g. Hence the non- leading 
part in terms of the function if) (z) is the crucial part that contributes in physical processes. 

A more complicate type of quantity that we needed in Eq. ()3.39|) takes the form of a double 
integral 

/ 

vk~ 1/2 F (it) k F (it) k~ iI2 i 



dn / dn 
ii Jo 



tanh ^ 



i cosh" 



%F(tanh 2 (f)) 



[K-n 



\ 



F(tanh 2 (^)) 



tanh ■ 



ftcosh 



\ 
/ 

(4.16) 
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where 

f/0__, = V ° ^ oV ° ( K "> = o < K l L oK > +9 < K l L o| -/«'>■ (4.17) 



The form of (k ) kk > in terms of matrix elements of Lo in the k basis is shown in the next 
section. We apply this formula to compute 7 (0) in Eq. ()3.39|) . The integrations in Eq. ()4.16| ) 
are well behaved for F (it) = (1 + it)~ 1 . But the infinite sum [k ) kk , is shown to diverge 
in Eqs. fl5.36M5.40J) . It is evaluated there more carefully by inserting a regulator in the form 
[e~ TKo k ) , . It is finite at finite r, but when r is small it behaves like 1/r times a quickly 
oscillating factor. The quick oscillations are not sufficient to overcome the divergence of the 
1/r factor. Therefore, this leads to the divergent result for 7 (0) in Eq. ()3.39|) which we used 
in the previous section. 

These integrals seem difficult, but can be easily performed with an algebraic computer 
program, such as Mapple or Mathematica. For example we obtain 







tanh — 2\/2 

dnVAK) ^-5 — J ^ = 0.544 33 =—-^ (4.18) 

; 3 + tanh 2 f Vk cosh 2 ?f 3^ 




, s 2 tanh ^ / tanh ^ 2 

cIkV 2 (k) i— J ^ = 0.222 22 =- 4.19 

3 + tanh 2 ^ V Kcosh 2 ^ 9 v ; 



' tanh^ 


2 


k cosh 2 ^ 3 - 


h tanh 2 ^ 



tanh^ 27 

^ = 0.523 25 = In — (4.20) 



k cosh 2 22* 16 



We now turn to the more specialized cases of relevance for our computations in this paper. 
We would like to evaluate the scalar u = VKo ^jt^^g > the vectors b Q = (^^f^o^vj 

and b e = {t^ii^o 1 ^ 2 vj and matrix elements M OQ i = (^§) oq/ , M ee i = (^f|) ee ,, M oa / = 

( 1+3& ) " ee ' = (i+is) '■ ^ a PPlyi n g the above integral formulas we obtained easily the 
following values which are the only ones we actually needed in this paper to perform the 
expansions of a (r) , (r) , 7 (r) , a 4 (r) in powers of e~ r up to order O (e~ 6r ) 

2y/2 , 2 , 22^2 , 1V2 , 2 3 / 2 67 
h = -^-f=, b 2 = -, h = 64 = h = 5—= 4.21 

3^ 9' 3 3 5 ' 4 3 5 ' (27)'^ 



Mu = — , M 22 = H M13 = -^5, a; = In — , (4.22) 



11 w 19 ~ 2 4 5 v / 3 

_ M 22 = — — -, Mi 3 = 

27' 9(27)' (27Y 



M n = — , M 22 = 777 ^ T , M13 = — (4.23) 



We also show how to compute some determinants by using Eq. (|4.14jl . We especially make 
a point to separate the contribution associated with the regulator iV which will be sent to 
infinity later. The steps of computation are shown for the following explicit example 

det ((«)) = e ^K2Ve*)/;°ln(tanh^)d K x ^ )+f(-f ))) ln(t^ ?f)d K 
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(2iVe A ) 1 (0.79788) = (2iVe A ) 3 J- = (vriVe A ) 



-1/2 



(4.25) 



As seen, the finite factor that comes from the non-leading term involving ijj is non-trivial. 
We will show later in Eq. (|5.fi2j) how both the leading and non- leading terms agree with the 
regulated discrete basis. In a similar way we obtain 




(2Ne A ) 18 (0.980 52) 
(2Ne A y^ (0.91746) 

(2Ne A y^ (0.993 40) 
(2Ne A )~^ (0.905 34) 



(4.26) 
(4.27) 

(4.28) 
(4.29) 



Note the magical number 26 that popped up as 13 = 26/2 in the last ghost determinant. 
Inserting these determinants in the expression for 04 in Eq. ()3.5|) we see that the N dependence 
drops out at r = 0, only for d = 26, 

(2Ne A )~^ (0.905 34)' 2 



a 4 (0) = g< 



27 
16 



9 



27 
16 



(2Ne A )~™ (0.993 40) 
(0.973 63). 
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(4.30) 



(4.31) 



So the result is independent of the cutoff N . A . This is an example of the magical role of 
d = 26. Note that a^, as given initially in [13j[l5j has factors of additional determinants 
that are actually divergent, but those get absorbed into the definition of g and produce the 
overall finite factor g 2 (27/16) 3 /4 (multiplied by another factor of 4 because of 4 diagrams). 

In the case of the bare coupling in Eq. (j2.3j) . inserting the results above we get 

3 / 2 / 1 \ -(d-6)/4 



90 = l 9 {T6) {( 2NeA y (°- 797 

= gc (2Ne A ) 5 for d = 26 , with c 



(2N)~* (0.980 52) 



27 x3/2 



16 



(2.203 0) 



(2N)-v (0.91746) 
(4.32) 



We see that the bare coupling diverges as (iVja. 

In the Appendix we have included many more values for the higher modes of the vectors b 
and matrices M, M. Although these were not needed for our analytic calculation, they were 
used to obtain the order ten result for F(x). Also the general formulas we give there are 
useful to compute to arbitrary powers of e~ T if one needs to refine our work in the future, 
or they could be useful in other applications. 
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V. MOYAL BASES AS THE j = REPRESENTATION OF SL(2,R) FOR L ,L ±1 

In this section we will study certain properties of the j = representation of SL(2, R) in 
relation to the continuous and discrete Moyal bases. The connection of SL(2, R) generated by 
Lq, L±\ and the continuous basis is known through the work in 21 1 1^. The relevance of the 
j = representation is emphasized in . Here we will study the j = representation with 
a different technique, mainly by focusing on the exponentiated group element e~ rL ° . We will 
extract some of the relevant properties of the j = representation that are needed in this 
paper. There are overlaps between our results and those in related papers [l^ [3| [3| 0] j^- 
We will see that all the computations in the previous sections, and similarly more general 
computations, amount to various matrix elements of functions of L , L±\ in the special j = 
representation of SL(2, R) . 

Let us first recall why we must focus on the j = representation. The Virasoro operators 
Lq, L±i are the generators of the SL(2, R) transformations on the open string basis X (r, a) . 
It is well known that this representation (in terms of differential operators) forms the zero 
Casimir representation of SL(2,i?) . Thus the oscillator basis labelled by the integers n = 
1, 2, 3, • ■ ■ is simply the case of the discrete series of SL(2, R) labelled as \j, m >, m — j + 1, 
j + 2, j + 3, ■ • • , for j = 0. 

We will use the Hermitian combinations Q\ = \ {L\ + L_i) and Q% — \ {L\ — that 
form the SL(2,i?) Lie algebra as follows [Q\,Q2\ = —iL , [L ,Qi] = iQ 2 , [L , Qq\ = —iQ\. 

Thus, the odd/even discrete bases in MSFT correspond to the eigenstates of the Virasoro 
operator Lo\n >= n\n >, n > 1, with n — o,e. The continuous basis is given by the eigenstate 
of the Virasoro operator Q\ = \ {L\ + satisfying Q\\k >= ||/« > . The change of basis 
is given by the functions < o\n >= 4gK ( K ) an d < e\n >= (/t) . The eigenvalue of Qi 
can have any sign, but in MSFT it is possible to make transformations to the range k > 
thanks to the symmetry property V Q (—k) = V Q (k) and V e (— k) = —V e (k) of these functions. 
We consider these as orthogonal matrices V OJ V e with matrix elements (V ) 0K = V (k) and 
(V e ) eK = V e (ft) , for k > 0, which satisfy orthogonality relations 

( V oVo) 0l02 = 5 0102 , (V e V e ) eie2 = 5 eie2 , (V V ) kk , = tfW (« - k') , (V e V e ) KK , = 6^ (k - k>) . 

(5.1) 

where a bar on a matrix means the transpose of the matrix. Here 5^°'^ (k — k') is basically 
the usual Dirac delta function, except for some extra care when both close to 

(see below), while their sum yields the usual delta function 5^ (k — k') + 5^ (k — k') = 
25 (k — k') . These are equivalent to the orthogonality and completeness relations of the 
discrete and continuous bases 

< n\n >= 5 nn >, < k\k >= 5 (k — k) , |n >< n\ = 1, / oIk\k >< k\ = 1. (5.2) 

n=l 

With the normalization given above the functions V Q (n) , V e (k) are identical to the functions 
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\[2v (k) , \plv e (k) in 21] 16], including an overall normalization of a/2 (consistent with 
counterpart orthogonal transformations in the discrete basis as explained in 0]). They are 
given by the following generating functions (in the range k > 0) 

EV (k) (tanz)° sinh (kz) ^-^ V e (k) (tan z) e cosh (kz) — 1 
T = / => y. 7= = i = I 5 - 3 ) 

v° J«;smh(^) e ve </« s inh(f) 



To define the functions for negative values of k the right hand side of these equations should 
be multiplied by the sign function e (k) for consistency with the symmetry properties of 
V (—k) = V (k) and V e (—k) = —V (k) . 

A fundamental quantity in MSFT is the matrix T eo [9|. In computations T usually 
appears in the form t eo as given in Eq. (j4.3|) . This matrix is diagonalized [3] by t eo = 
J °° dnV e (k) t (/t) V Q (k) with the eigenvalues r (k) = tanh (ttk/4) . To understand this matrix 
we introduce the notion of the operator t. Then t eo can be recognized to be just the matrix 
elements of the operator t in the discrete basis 

ttQ 1 1 — e^ nQl 
t = tanh = t eo =< eltlo > . (5.4) 

This is proven by introducing identity in the k basis and writing 

POO 

t eo =< e\t\o >= / dn < e\t\K >< k\o > (5.5) 



/°° 1 UK 1 
dn—=V (k) tanh —V e (k) (5.6) 



7TK 







cZkT4 (ac) tanh («) , (5.7) 



where we have used the symmetry of the functions in going from the second line to the third. 
By the same argument we can show t ee > = and t QO i = 0. So the matrix representation of 
the operator t is block-off-diagonal in the even/odd basis. 

Now we note the special role of the k = state and its relation to v Q . The vector v Q is a 
crucial state in MSFT as seen in the calculations in this paper. Its special role in relation 
to anomalies first emerged in [l^] as a state closely related to the midpoint of the string. 
From the generating functions in Eq. (|5.3j) it is seen that at k = the even functions vanish 
V e (0) = 0, while the odd ones V Q (0) are finite, 

K(0) = 0, V o (0) = l-V^b~v o . (5.8) 



where v Q is given in Eq. ()4.3|) . It is a normalized vector ^ 0>1 v o v = 1 that is a zero mode of 
the matrix T : ^ 0>1 T eo v = 0. Equivalently, V a (0) is the zero mode of the matrix t eo . This 
is easily seen from the expression Eq. (j5.4j) for the t operator acting on the zero eigenstate in 
the kappa basis Qi|0 >= 

t eo V;(0) = < — > e- nQl \0 >= |0 > . (5.9) 
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Note that we can write v a in terms of the k = state as follows 



= (2/v^F) <o\ — 10> . (5.10) 
Lq 



We can use the above interpretation of the operator t and and its zero mode v Q to give an 
SL(2, R) operator representation of the computations we performed in the previous section. 
In particular, we see that b D ,b e of Eqs. (|3.8|) can now be written as matrix elements in the 
j = representation of SL(2, R) 

2 2 1 2 2t 1 

b = -= < o\— — 1T - 10 >, b e = -= < e\— — 10 >, (5.11) 
V 71 S + t z L V 71 " 3 + ' 

Note that < e|TTi2T-|0 >= since < e|0 >= 0. Hence we can remove the bras and write the 
states 

^ >=^3TF> >• |6 < >= ^sTF> > (5 ' 12) 

where now the o, e indices may be interpreted as an odd or even number of powers of the 
operator t. Similarly, the quantities a (r) , (3 (r) , 7 (r) that appear in the 4-tachyon amplitude 
in Eqs. fl3.2H3.4|) become expectation values in the k = states. For example 

f3 ( r ) = bP {o) e~ TK (1 - Me- Tft Me~ TK ) ^ P {o) b, (5.13) 

= 1 < 01— — e- rLo ( 1 - 1 ~ * e- TLo1 ~ * 2 e- TLo ) ' — -10 > (5 14) 

Once we notice these forms we are tempted to use the commutation rules of the operators of 
SL(2, R) to simplify and compute these quantities. In particular note that e~ w< ^ 1 L e^ 1 = 
Lo cos a; + 1Q2 smu gives the following interesting properties of t and e -71 "^ 1 for u = tt 

e~ wQl L = -Loe-^ 1 , tL = L r\ —t = r 1 -^, t— = — r 1 , etc. (5.15) 

Lo Lo Lo Lo 

Using these it appears as if we may perform a number of simplifications very efficiently Alas, 
there is the problem of anomalies precisely because of the existence of the k = state on 
which t vanishes. So its inverse cannot be used in a cavalier fashion. In fact not all of the 
above formulas are necessarily valid in general since they depend on which state they are 
applied. When we introduce a complete set of intermediate states between operators, the 
k = state must appear, causing problems. One can easily generate inconsistent results 
from the same formula with naive manipulations of t. Examples of inconsistencies can be 
easily constructed by taking matrix elements and using naively the formulas above on the 
k = state. The type of formal manipulations suggested above need to be justified at every 
step or otherwise replaced by the correct result. 

Thus, practical computations in string field theory require a regulator at k = because 
the matrix t eo has certain associativity anomalies that come from the zero modes reviewed 
above. This is precisely the zero mode of T eo and the associativity anomaly issue which are 
intimately related to midpoint issues in string field theory as first understood in |IoT |. 
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A regulated version of the matrix t eo with many nice mathematical properties was in- 
troduced in [3 0] • The regulated matrices retain the properties of the operator t that 
are always valid (see Eq.(2.15) in 111]). In the regularized theory, the continuous func- 
tions V Q (ft) , V e (k) can be thought of as the large N limit of N x N matrices (V ) ok , (V e ) ek 
with discrete values of k which label the N eigenvalues of the regularized N x N matrix 
t eo = V e rV . The large N limits of the matrix elements {V ) ok , (V e ) ek become the functions 
V (k) , V e (k) , and the discrete eigenvalues become — > r (k) = tanh (ttk/4) . In the limit 
iV — > oo, summation over discrete k is replaced by J °° and the discrete delta function Skk' 
is replaced by the continuous delta functions 5^° ,e)> (k — k') with non-negative k, k! > 0, but 
with some extra care when both close to zero, as we will see below. 

In this section we try to relate the large TV limit of the regulated theory in the discrete 
basis to the regulated continuous kappa basis, and from this extract computational methods 
that take advantage correctly of both the discrete and continuous bases. For this purpose 
we will study several quantities by taking advantage of the SL(2, R) group generated by the 
Virasoro operators L , L±, L_i 



A. Propagator < n\e 7 °\k > 

The quantity D KK i (r) =< K,\e~ TLo \n' > is a representation of a group element of SL(2, R) 
with Lie generators (L , L ±1 = QiT 1Q2) > since |« > is the basis that diagonalizes the gen- 
erator Qi — ^ (Li + L-i) , with Qi\k >= > . Its integral and derivatives are quantities 
of interest in various explicit computations in MSFT 

1 f°° 

< k\ — \k' >= / drD KK , (r) , < k\L \k' >= [-d T D KK , (r)] T=0 . (5.16) 
L o Jo 

In this section we compute D KK i (r) as a group element of SL(2, R) in the special represen- 
tation j = 0. Of course, in the discrete basis e~ rL ° is diagonal, therefore we expect 



D KK , (r) =< K\e~ TL ° I \e >< e\ + ^ \o >< o\ ) \k > (5.17) 

\ e o 

= \ (V e-^V ) KK , + \ (V e e-^V e ) KK , . (5.18) 

where the diagonal matrices K e = diag (2, 4, 6, • • • ) and k q = diag (1, 3, 5, • • • ) , and 

the V Q , V e are treated as matrices with matrix elements that are the functions (V ) OK — V (k) , 
etc. 



Following the methods in [22j, Chapter 7.4 (leading to Eq.(7.65)), we can derive a differ- 
ential equation for D KK > (r) as follows. We use the fact that the eigenvalue of the quadratic 
Casimir operator vanishes for the states \k > to write 

< K\e- rLo ((Lq - Ql - Ql)) \k >= 0. (5.19) 
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The first two terms can be rewritten as follows 

< K\e- TLo L 2 \K' >= d 2 D KK , (r) , < K | e - TL ° (-Q 2 ) \k' >= - D KK , (r) . (5.20) 

The third term < K>\e~ TLo (— Q|) \n' > can also be computed by noting the following property, 
e TLo Qie~ rL ° = Qi coshr + iQ 2 sinhr, which can be rewritten as 

e-^HQ, = -J—Q^o - ^L e ^ Q (5 21) 
sinn r sinh r 

Multiplying Eq. ()5.21|) with iQ 2 , using Q1Q2 = Q2Q1 — iL , and using Eq. (j5.21|) again, gives 



r - TLo Ql = (e- Ti °Q 2 ) - ( e ~ TL °Qi) Qi - ^^Lo (5.22) 

smhr sinhr sinhr 

1 ^1 n-r coshr ^ * ^ cosh 2 r , ^ 9 coshr^ . 
-Qle- TLo - 2——^—Q 1 e °Q 1 + -— 2~ e <2? + ; ~ d T & 



sinh r sinh r sinh r smhr 

(5.23) 

The matrix elements of this relation are easily evaluated in terms of D KK i (r) . Combining 
the three terms in Eq. (j5.19|) we obtain the differential equation satisfied by D KK i (r) 

o2 coshr k 2 + (k 1 ) 2 - 2kk' coshr . _ 

2 + —r—d T + ; , -D KK ' t = 0. 5.24 

smh r 4 srnlr r ' 

By taking D KK , (r) = 2 (1 - z)"*^"')/ 4 F (z) , with 



sinh -J = -4e" T (l - e~ T ) , (5.25) 
this differential equation becomes the hypergeometric differential equation for F (z) 

d \ 2 c— (l + a + b)zd ah .,, 

m) + la-,, fc'^) |F(;) = - (5 ' 26) 

with F (a, 6; c; 2) ~ hypergeom ([a, 6] , [c] , z) , and 

a = l-z|, 6 = 1 — i^-, c = 2. (5.27) 
We need the solution that satisfies the boundary conditions that follow from Eq. (j5.18|) 

D kk ,{t) -> ^Vi («) VI («') e" r , D KK i (t) -> ( J( /c -«') j (5.28) 



where Vj («) = ^ sin £™ > as obtained from Eq. (j5.3|) . Therefore, the desired solution is given 
by 



, l-e~ r \ 2 2 / m m' -4e" T 
A<*/(r = J—— —7— -o F l 2: 



IK I IK 



"'" ' ' l/ sinh f sinh ^2(1 -e-) 2 VI + e"V \ 2' 2' ' (1 _ e -) 2 

(5.29) 
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1 — zf a b F (c — a, c — b; c; z) , we see that 



By using the linear identity F (a, b; c; z) - 
the expression for D KK > (r) is real. At r = oo the boundary condition is satisfied since 
F (a, b; c; 0) = 1. For small r ~ 0, and any k, K f , we use the asymptotic behavior of the 
hypergeometric function to write 



1 



KK 



Dkk ' ^ t~o 2 3 V sinh 2£ S inh ^ 



r i + 



ft') In 



r(i + f)r(i-f) 



+ c.c. 



[K — K' 



(5.30) 



where c.c. is the complex conjugate. For small (V — k) ~ and small r ~ we can write 
r (l + i§) r (l — l f) = ?f (sinh ^f)" 1 , so that the singular behavior of D KK > (r) , consistent 
with the delta function, becomes evident as follows 



sin ( (k — k') In 



7T (« — K 1 ) 



2tt 



in 



in 



We also record here the special values for k, k' = 0, that follow from F (a, l;2;z) 

{l-z)- a z-(\-z)- a +l 
z(l-a) 



at any r 



< K\e 



-tLq 



< 0|e- rLo |0 >= D 00 (t) = — In 



sin 



>= D 



2 I 1-e" 



^2™ sinh ?f 



1 



1 + e" T 



2tt V 1 - e_ 

where the second line follows by setting k = in the first line. 
The integrals of these quantities give < «|^|0 > and < 0|-^|0 > 



< K\ — 10 > = 



drD K0 (t 



7r tanh ?f 
16k cosh 2 £f 



1/2 



, 1 . 7T 



This agrees with an independent calculation in Eq. (|5.fi4j) . 

The derivative near r ~ gives information about matrix elements of L 

< k\L \k' >= [—d T D KK i (r)] T _ 

"r(i + 



KK 



2 3 A/ sinh sinh ^ 



{(k—k') \ re')ln 
2 1 e 



r(i + f)r(i-f)r 



+ c.c. 



(5.32) 
(5.33) 



(5.34) 



(5.35) 
(5.36) 



This diverges due to the factor 1/r. More simply, the quantity < k|L |0 > is obtained also 
by studying the derivative of Eq. (j5.32J) at any r 



< K\e- TLo L \0 >= -d T D 



Ke 



-t COS 



-ln , 

2 111 \ l-e- 



n0 



1-e 



-It 



^JIixk sinh ^ 



(5.37) 
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Near r ~ one sees that < «;|L |0 > and < 0|L |0 > have a singular behavior 

f m \V 2 cos [fin , , 1 

< ,|L |0 >^ (-^) L | _UJ ) < 0|Lo|0 ^ _ (5 . 38) 

Let us also compute < k\Lq\k' > for k' ~ K from Eq. (|5.3(jf) or from Eq. (J5.31j) . We obtain 
the singular behavior 



cos 



[< «I^o|k' >U„ = [-flr^ (r)]_ ,^ K « (5.39) 

which is consistent with Eq. (|5.38j l. but now gives 

< >= [ d T D KK (r)] TM) « (5.40) 

Z7TT 

at any «. The divergence of < k|L |ac' > is also discussed in jig . 



The quickly oscillating factor exp(i (k — k') In y 2) j n Eq. (j5.36j) suggests that only the 
neighborhood of k' ~ k is relevant in a steepest descent approximation when the func- 
tion < «;|Lq|k' > is integrated with smooth functions J dn J dn' f («) < «|Lo|/c' > g («') 

~ 2^7 / ( K ) 9 ( K ) / dq cos ^gln y^)> where we have used Eq. (j5.39|) . and implied a cutoff 
in the integration in q which depends on the functions /, g. This expression is also divergent 
when r vanishes for well behaved behavior of f,g. This last form relates directly to the 
quantity 7 (0) in Eq. ()3.39|) through Eqs. (|4.1fi|4.17|) as the regulator r is removed. We used 
the singular behavior 7 (0) = 00 to determine the mid-range slope (i.e. differentiability) of 
the off-shell factor in the 4-tachyon scattering. 



1. Computation of (V e e KeT V e ) KK , and (V Q e K ° T V ) , 

To compute (y e e~ KeT V e ) , and (V e~ KoT Vo) kk , , we first recall their sum as in Eq. (j5.18j) . 
and then extract the individual terms from the symmetry under reflections with respect to 



k or k'. 



(V e^ T V ) KK , = D KK i (r) + D { _ k)r , (r) = D KK , (r) + D« ( _ K ,) (r) , (5.41) 
(V e e-^ T V e ) KK , = D KK , (r) - D ( _ kV (r) = D KK > (r) - D re( _ re0 (r) . (5.42) 

where we insert our result ()5.29|) for D KK i (r) . 

In particular, if we take k'—0 we get (V e e~ KeT V e ) kQ = 0, while (K>e~ K ° r V o ) K0 = 2D k0 (r) , 
which implies 

v K X J 2 e-^V (k) = 4= (Ke- K ° r K) K0 = 4= < 0|e- Lor |«: >= -^D K o (r) (5.43) 
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Taking also k = gives 



v K e-*° T v = - {V e~*° T V ) m = -< 0\e- L ° T \0 >= -D m (r) (5.45) 

We are now ready to tackle the quantity v k l ^ 2 V (k) = -7= J2o>o ^> (0) o^> ( K ) wrnc h 
is needed in explicit computations in MSFT, as in Eq. (j4.7|) and the equations that follow 
it. Since V e (0) = 0, only the odd V Q (0) contribute when the identity 1 = J2 n \ n n \ ls 
inserted in < k|t-1|0 > 

1/1 \ , 1 , / 7r tanh 



V -V =< k\ — \0 >= f— , (5.47) 

2 V °o Ao V ^cosh 2 ^ ' 

where we used the result in Eq. (|5.H4jl . From this we obtain immediately 

This result is used to relate regulated discrete basis computations of the type 



V n K, 



1 2 {f (tt)) oo , k o ^ 2 v i in the large N limit to kappa basis computations as shown in 



the previous section. 



2. Computation of p(n) =< k\k > 



We have argued in Eq. ([5.31|) that near r = and k' ~ k we obtain D KK r (0) =< >= 
5 (k — k') . The quantity p (k) =< k\k > appears in the evaluation of determinants as in 



det (f(tt))= exp 



/ dnp (k) In / ^tanh 2 — ^ 

J 



(5.49) 



Since we expect a singular behavior, p (k) ~ 5 (0) , we would like to perform such computa- 
tions with a finite number of modes 2N and send iV to infinity at the end of the computation. 
Thus, on the left side we have a finite determinant that depends on N, which in principle 
can be computed in the regulated basis, at least numerically. To perform such careful com- 
putations analytically by taking advantage of the k basis, we define p («) =< k\k > carefully 
at k' — k for small but finite r, and compare to a similar calculation with the cutoff in terms 
of the number of modes 2N. Thus, for k' = ±k the cases D KK (t) , D K (- K \ (r) , as r —>■ 0, are 
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used to compute the quantity p (k, t) ~< k|k > as follows. From Eq. (J5.3U)) and from the 
asymptotic behavior of the hypergeometric function at large (— z) 

F (a, a, c; z) - F '^7^ , (In (-z) + 2 t/j (1) - ^> (a) - ^ (c - a)) , (5.50) 
1 (a) 1 (c — a) 

we obtain 

A. (r) - = ^ (in ( I ) 2 + 2 (1) - (l - I ) - (l + I ) ) + O (r 2 ) . (5.51) 

D. ( _ w) (r) -> = , L _ (! + »«) ■ <wln i _£0^) e -«m|V 0(r2) (552) 

K(K)v; r~o 2 3 ismhf ^r 2 (i + f ) r 2 (i-f) y v ; v ; 

Note the tricky behavior near k = : as long as k is not close to zero, the quantity 
lim T ^ (r) oscillates rapidly and becomes negligible as a distribution compared to 

lim T ^ D KK (r) . However, near k = we get 

1 /2\ 2 / vrm 2 - In 32 - \ 
B <-* < T > ^~.= S [r) + K ( "-24 1 -1ST + a ° 4783 + ° ^ (5 ' 53) 



Hence, if k 7^ not small, and r is small, we neglect D^ K ) K (r) and get (V e TO K) 
(V e e -Te y e ) ~ D Kfi (r) ~ p («, r) with 



KK 



1/. /2\ 2 , /N , / ifc\ , / m 



p {k,t) = — lnl-l +2 ^(l)-^M--j-^M+-) ). (5.54) 



But, close to k = 0, because of Eq. (j5.53j) . the behavior is 

k,t~0: (K ) e-™K ) ) 00 = 2[p(0,r)U = |-ln^ , {V e e~^V e ) QQ = 0. 



(5.55) 



The factor of 2 computed here is important to establish consistency between the regulators 
in the discrete and continuous bases. 

We concentrate on p (k) for the general k. We already understand its dependence on k 
explicitly, and therefore write it in the form 

P («) = P (0) + ^ ( 2 V> (1) - V> (l - J ) - V> (l + J ) ) • (5-56) 
We know p (0) is divergent, and is given by 

P (0) = I (VoV o ) Q0 =< 0\0>=l (v KoVo) regulated = i (In (2iV) + A') (5.57) 

where we used a finite number of modes 2N as the regulator and the constant A' depends 
on further details of the regulator. Therefore the full regulated p (k) is given by 

p («) = i- (In (2iV) + A) - i- ^ (l - y ) + V (l + y ) ) (5-58) 
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where A is a constant. Note that this constant can be absorbed into a redefinition of N 
by writing In (2iV) + A = In (2iVe A ). We will further determine A = by comparing to a 
particular regulator in the discrete basis of MSFT in which we keep the frequencies as simple 
as possible, namely K n = n. In principle it is a dangerous business to compare diverging 
quantities when using different regulators. However, as we have already mentioned, since 
the N dependence will drop out in certain computations in d = 26 then A will also drop 
out, and hence its value may not be important. 



B. Exact analytic tests 



As tests of the formulas we derived above in the kappa basis we compare to exact compu- 
tations in the regulated 6 discrete basis, in the limit iV — ► oo. The first case is vk q 1 ^ 2 1k 1 ^ 2 v 
where we compare Eq. ([4.13|) versus using directly the v D in 14.31 or the regulated one in 
footnote E3 The two ways of computing give the same answer: 

K- 1/2 i«» I/2 «)*«« t « = £ E j^rf + « lniV > /*") ^ ^ = °- 85256 < 5 - 59 > 

roc tanh — 

(vk~ 1/2 1k' 1/2 v)^ = dn 1 = 0.852 56. (5.60) 

V ° ° >kappa Jo K (coshf) 2 1 ' 



The second case is det(tt) = exp(/ °° dnp(K)\n (f (tanh 2 ^))) where we use p (k) in 
Eq. (j4.24|) . while the exact computation in the regulated discrete basis was given in [lG| 
for any N as 



11 



. V 2 detn l-3-5---(2iV-l) 



det {tt) discreu = (1 + ww); egulated = — ^ 2 . 4 . 6 ... (2jV) (5.61) 
" T(2N + 1) ~{«N)-*(l + ±Y'\ (5.62) 



2 2N T 2 (N + 1) v ' \ 2N t 

where the last expression is a good approximation for any N > 1 (including small N). Note 
both the leading N behavior as well as the overall factor (n/2)~ 1 ^ 2 (2N)~ 1 ^ 2 , both of which 

6 To avoid anomalies we use the standard regularization of MSFT If we make the choice K n = n, 

at finite N the regulated w e , v a , T eo , R oe take the form T eo = (f?/f?) , R oe = R ( o? ] (felfo) » 

w e = We fe i v o = v o°°^ I fo where T^°°\ R^°°\ u/°°) , are identical to the expressions in the infinite 
limit as in Ea. H4.3|l but truncated as N x N matrices. The factors /„ for n = o, e are given by /„ = 

J r^+i-C)r(jv+i^^ ' r( ( N+i ) ) " Tllis was com P uted b y simply inserting « n = n in the general regulated 
formulas for arbitrary n n given in ^(J . Note that the deformation factor f£ stays pretty close to 1 

(in the range 1.00-1.075) when 1 < n < N, even for finite N, but grows as n approaches 27V. For n = 2N, 
and large TV, we get f^ N ~ (^j^) 1 ^ 4 ■ The strong deformation at the edges of the matrix can be avoided 
by using a different function n n as a function of n. However, in our experience the simple choice n n = n 
seems to work well in numerical computations even at small values of N. In the present example we see 
that it also works exactly at infinite N in comparison to the continuous basis. 
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are significant for the comparison of regulators in different bases. We compare this to the 
computation of det(tt) fe a in the kappa basis using p(k) as given in Eq. (j4.24|) . We see 
that the leading term ^- In (2iVe A ) in p (k) reproduces the correct N dependence, while the 

second kappa dependent term produces precisely a factor of \ leading to the total result 



det (it) kappa = (vrA^e A ) ^ 2 . Comparing to the regulated discrete basis result in Eq. (j5.62j ) 
we learn that we need A = to get agreement between the two regulated results. Then, we 
seem to have an agreement between the computational procedures in the regulated discrete 
basis and the regulated kappa basis. However, we must warn the reader that A = is still 
tentative because it is possible to change the regulator in the discrete basis, and we have not 
understood yet the principle that could fix it in either the discrete or the continuous basis. 
It is however significant that we have seen examples where N as well as A cancel together 
in finite quantities at the critical dimension d = 26. So, perhaps the value of A is not crucial 
as long as one is consistent in using the same regulator everywhere. 

The third case is the exact computation of vv in the regulated discrete basis for any N 



f- \ f ww\ det k 

W— . = (TT^A^ = 1 -d^ (5 - 63) 
k7=u 1 ~ 2 2Jv r 2 (n + 1) ^v-oo L ( 5 ' 64 ) 

At iV = oo in the kappa basis we use V (0) = \^/twv and obtain agreement with the 
discrete basis as follows 

= = ^ K >= 1 (5 ' 65) 

where we used Eq. fl5.34j) . 

In general we have noticed no problem in agreeing between the two bases for finite quan- 
tities where the effects of iV disappear. For such quantities it turns out that even a few 
modes in the discrete basis (N ~ 5—10) gives very reliable numerical estimates of the exact 
N = oo values. However we advise care when the quantity blows up or vanishes as a power TV 
(such as determinants). Finite iV estimates are not predictably good for such quantities and 
furthermore they seem to be cutoff dependent (such as A). When divergent quantities occur 
in combinations in which iV cancels, the various computational approaches seem to become 
reliable again, and furthermore in our experience, the result seems to be well estimated again 
at relatively small values of N. 

VI. NONPERTURBATIVE LANDSCAPE 



The perturbative computations above have implications for nonperturbative computations 
in MSFT such as the ones in jl^ . There a proposal was made for computing the true vacuum 
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of Witten's string field theory at the classical level. This involved solving the equation of 
motion that follows from the action in Eq. ljl.lj) 

(L - l)A + a'g A*A = , (6.1) 

for fields that are independent of the midpoint (hence the D 2 5 brane vacuum). This could 
be done by treating the energy of the midpoint (called 7) as a perturbation and writing 
the rest of L as a special star product form £ * A + A * C . All solutions of Eq. (j6.1|) . 
including the vacuum, were obtained exactly in the absence of 7. Then it was possible to 
setup a perturbation series in the midpoint energy 7 and compute analytically each term 
order by order. The lowest term was obtained as an exact solution A^ (x e ,p e ), where 
A(°> (x e ,p e ) turned out to be related, up to an energy dependent factor A^ = — (2/a'g) £0* 
P, to the twisted butterfly projector P = A^ (x e ,p e ) written in the MSFT basis 7 . The next 
perturbation A^ was also computed explicitly, while higher orders A^ n > could be computed 
in a straightforward way with similar methods. Finally the energy of the vacuum and the 
tension of the D25 brane were given analytically up to second order in 7 

T 25 = ±-S (A® + eA® + e 2 A^ + ■■■) (6.2) 

^25 

~^ + 4^+0(^ 2 ), (6-3) 



with 



«% 2 V3 2 

^^-^fe^-E^V *=(d-2)%^. (6.4) 
2 4 \ / — i ^-^ I 1 + ww 

\e>0 o>0 / 

When this result was obtained the divergent nature of v and 5 were confusing. However 
the divergence of the bare coupling constant g was not noticed. We have seen in this paper 
in Eqs. fj2.3l4.32p that go = gc (2iVe A ) 2 . A quick look at is, 5 show that they are both linearly 
divergent with (2N) as can be seen in a naive level truncation up to (2N) by inserting k q = o, 
n e = e and using the unregulated w e = —\/2(—l) e ^ 2 . Then the result for T 25 is perfectly 
finite at each order of the midpoint energy 7 since the factors of (2iV) 3 cancel between 
numerator and denominator in Eq. (j6.3|) . Thus, this nonperturbative computation should 
proceed by using the same regulator consistently, keeping only the leading terms to order 
(2iV) 3 in the numerator and dropping everything non- leading. It reminds one of the large 
N matrix computations. Progress on determining the vacuum state and energy using this 
approach will be reported in a future paper. 

The observation above resolves another puzzle as follows. We recall the expression for L$ 
that appears in the action of MSFT in Eq. (|l.l|) . We display the version in ^| including 

ghosts L = L natter + hf ost , 

e>0 v e re 7 



7 A proposal based on the butterfly that has parallels to our program appeared recently 
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+ — ^2.^-1^(2^^) -2Z. K «' (6 - 5) 

e>0 e \e>0 / n=l 

\ e / \ e' e / n=l 

where 3 = — iy/2a'-^= = —i\j2a'-J^— represents the center of mass momentum. We had 

OX OXcm 

commented before in previous papers that the terms involving (1 + ww)~ are tricky. Even 
though they appear to vanish as iV — > oo they actually contribute a finite term because 
infinite sums cancel the zero. The midpoint energy 7 mentioned above is in fact just the 
piece of L proportional to (1 + ww) -1 . 

But how about the divergent term in the first line | (1 + ww) 0$ involving the center 
of mass momentum /3q? This divergence caused concern for some colleagues. Actually we 
can explain that this is the correct behavior of this term because otherwise there will be 
no center of mass momentum dependence in the large N limit in nonperturbative physics. 
Let us start with the cubic term g$Tr (A 3 ) and replace A with A = g^A' so that it takes 
the form g^ 2 Tr (A') 3 . Doing the same to the quadratic term, we rewrite the action with 
an overall 2 as S = -g 2 (Tr ((l/2a') A' (L - 1) A' + (1/3) Tr (A') 3 )) . Now it has the 
form of the computed nonperturbative energy in Eq. (j6.3|) and, from the discussion we gave 
above, we see that the numerator must behave like iV 3 to compensate for the behavior of 
go — eg (2iVe A ) 3//2 in the denominator. It is clear that each A' accounts for a factor of iV in 
the numerator and also that L contributes a factor of N in the quadratic term. Now, if we 
also include momentum dependence (or x dependence) in L Q the only way that momentum 
will not be negligible is by getting the help from the factor (1 + ww) ~ 2N in the form that 
it appears in Lq = • • • + 1 (1 + ww) (3$. This explains that this factor has precisely the correct 
behavior, and how it contributes in nonperturbative phenomena. 

Of course, there is a way of eliminating the confusing factor by renormalizing A 1 = 
(l + ww)A and defining L = (1 + ww)~ x L , and further absorbing the extra factors in 
a definition g = (1 + ww) 3 ^ 2 g . Then in the new L every sum is divided by the factor 
(1 + ww) ~ 2N which is reminiscent of the finite quantities Tr (M) /N in large iV matrix 
theories. If the theory is redefined in this manner in terms of (A, L j then every term of 
the action in every computation should be finite at every step, just like the leading terms in 
large N matrix theories. 
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VII. CONCLUSION 

In this paper we have demonstrated that MSFT not only agreed in great detail with other 
computational approaches to Witten's string field theory, but it also led to new results that 
were not obtained before. We have developed several practical and theoretical tools on the 
way to the new results, and we have indicated how certain nonperturbative computations 
can be conducted by using the information provided in the present paper. 
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IX. APPENDIX 

Anticipating future applications we compute also the following quantities by performing 
the integrals, which can be done easily by using an algebraic program. We give the results 
for five significant figures 



/ 22.222 \ 
-11.058 




(9.1) 



4.1275 



V 



/ 



/ 54.433 \ 
-12.804 




(9.2) 



3. 243 2 



V 



/ 
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(it) oo/ = 10" 



40. 528 
23. 399 

-14. 097 
9.8718 

-7.4881 



23. 399 
85.56 
10.464 



-14.097 
10.464 
91.684 



-8.136 5 6.850 6 
6.6112 -5.795 



9.8718 
-8. 136 5 
6.850 6 
94. 133 
5. 1061 



-7.4881 
6.6112 
-5. 795 
5. 1061 
95.46 



V 



\ 



(9.3) 



/ 



Note the increasing diagonal, although the off-diagonals are much smaller. The increase on 
the diagonal (it) 00 quickly slows down as o increases, and stays around 160 < (it) 00 < 167 in 
the range 100 < o < 175 and then makes a sharp drop reaching (it) 00 —2. 261 4 at o = 185, 
and continues to drop slowly as o increases. 

The following are exact computations, like the ones above, obtained directly by using the 
integrals, not by inserting a truncated form of tt or it. If compared to what follows from 
the truncated ti or it one finds results that disagree although they are in a similar range of 
values. Therefore, the contributions from the higher modes are not always negligible. 



tt 



3 + tt 



M ee , 



10" 



x 



/ 5.349 8 -3.678 7 2.753 3 -2.183 7 1.8013 
-3.678 7 2.8931 -2.333 7 1.942 2 -1.6571 
2.7533 -2.3337 1.9739 -1.6978 1.4843 
-2.183 7 1.942 2 -1.6978 1.4961 -1.332 3 
1.8013 -1.6571 1.4843 -1.3323 1.2039 

V ; ; ; ; ; 



\ 



/ 

(9.4) 



tt 



3 + tt 



M 00 , 



10" 



x 



/ 18.519 -7.603 4.725 9 -3.393 3 2.628 7 

-7.603 4.536 9 -3.288 2 2.572 3 -2.106 

4.725 9 -3.288 2 2.582 -2.124 3 1.8012 

-3.393 3 2.572 3 -2.124 3 1.808 8 -1.572 5 

2.628 7 -2.106 1.8012 -1.572 5 1.393 3 

V ; : : : : 



\ 



/ 

(9.5) 



tt 



\ + 3tt 



M PJ 



10~ 2 x 



/ 7.818 9 -5.748 4.536 7 -3.7613 3.223 3 
-5.748 4.657 7 -3.872 5 3.315 4 -2.903 5 
4.536 7 -3.872 5 3.326 6 -2.912 2 2.5918 
-3.7613 3.315 4 -2.912 2 2.5913 -2.334 8 
3.223 3 -2.903 5 2.5918 -2.334 8 2.1241 

V ; : : : : 



(9.6) 
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tt- 1 

1 + 3tt 



M 00 , 



10" 



/ 40.741 19.007 
19.007 10.664 
-12.905 3.288 2 



-12.905 9.9171 -8.1191 
7.800 8 -6.2376 5.233 4 
5.9573 4.895 7 -4.1851 



9.9171 -6.2376 4.895 7 4.1014 3.5561 
-8.1191 5.233 4 -4.1851 3.5561 3.1171 



Other matrix elements of interest are obtained from the ones above, such as 
M, and 2g±j£ = 1 — M, etc. 



/ 

(9.7) 

= 1 + 



We can also evaluate some of the quantities by using other methods and compare to the 
results above. In particular, by exact summation over the infinite modes we can express 
the infinite matrices {tt) ee , and (it) 00 , in terms of the generalized hypergeometric function as 
follows 



32 (i) c 



'OO' 



" tt 2 (o 2 -4) ((o') 2 -4) 

(\ o 
x hypergeom I 2, 1 + — , 1 



(9.1 



o' o' o 

— , 1 + — , 1 

2' 2' 2 



o 

2 + -,2 
2' 



o d 
-.2 ,2 

2' 2' 



2" 



The values of the hypergeometric function exactly agree with the results of the integrals 
given above. 

Also we recognize the quantities u, b a , b e , M 00 /, M ee >, M 00 >, M ee < from our earlier work on the 



computation of the Neumann matrices by using the Moyal product [11[ |15| and comparing 
to Neumann coefficients which were obtained from conformal field theory [7[. Therefore, in 
some sense we already knew the result for these quantities. However, what we knew is only 

and serves to confirm the 



a special case of the more general formulas given in Eq. ( 
general method. For example (it) 00 , (tt) ee , or the determinants evaluated in the text cannot 
be obtained by using just the Neumann matrices. Thus, from our earlier work we extract 
the following results given in terms of the vectors A e ,A ,B e , B whose numerical values are 
given by the following generating functions 



l + z 



1/3 



1 + iz 
1 



IZ 



2/3 



In terms of these we have the desired quantities as obtained from |TJJ [15] with a little algebra 



b 

be 



3 + tt 
2 



2 A ~ 
3^' K 



t ^K-^V 

3 + tt ° 



■y/2 



A r , 



t- 



3tt 
2 



1 + 3U 



2_B< L 
3^~o 

3 



(9.10) 
(9.11) 
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and 



1 - tt 



\{l-Al-2 J27=i K» + 2 E:"= 2 Al) , foro = d 



1 _ tf\ f 1(1 + A l + 2 C»=2 ^ - 2 E^=i ^») , for e = e> 



M 00 , 



"-Ijm^-i (913) 

r = ft_U = ( If^fL + (^ + ^U (9 14) 

\3+Uj„ i^Ti{ e + o + e-o ) [JAi> 

where a bar on a vector, such as A, means its transpose. Similarly we have 

1 _ tt \ \ -U 1 - A l - 2 E°o-i Al„ + 2 A 2 ,, - 2A oJ B c ) , foro = d 
— I = \ i / — ( (ab+ba) , (ab-ba) A 

- + 3«A , | -Iv^ 7 0+0 > - 1 _/"" j , foro^ d 

(9.15) 

1 _ tt \ ( -H 1 + A l + 2 C»=2 ^ - 2 E^li ^» - 2A e 5 e ) , /or e = e' 
= I = \ 1 / — / (ab+ba) , (ab-ba) A 

- + 3« Ae' \ ^ e+e >' - 1 e _> ] , for e ^ e' 

(9.16) 

y _ A 1 ^ _V^( (^-BA) eo (AB + BA) \ 
eo \l + m) eo AVS\ e + o e-o J' K ' J 

The generating functions give the following values for A e , A a , B e , B which are useful in the 
present paper 



t p 2 7 2 x 19 4 2 x 409 fi 2 x 11 x 283 8 2 x 220 721 in 

> A e z e = — -z 2 H — z A — z e H z 8 z w H 

/ v e 32 3 5 3 8 3 10 3 14 

(9.18) 

, n 2 2 x 11 , 2 x 67 , 2 x 1409 , 2 x 94993 q 

> A z° = -z — z 3 H —z 5 z 7 H z 9 H (9.19) 

3 34 36 39 313 v ) 

^ „ „ 2 3 2 2 4 11 4 2 3 523 „ 2 5 29 x 37 „ 2 3 323 381 in 

E = -3^ + -$r z - — z + 3T0 - -^14— - 10 + ■ ■ • ( 9 -20) 

„ 2 2 2 2 17 3 2 2 127 5 2 2 11 x 277 7 2 2 23 x 9839 9 . 

> SoZ = — 2 — z 3 H — z 5 - z 7 H z 9 H (9.21) 

3 34 36 39 313 v ) 

By applying these formulas we verified that the expressions in terms of the vectors A n , B n 
given in Eqs. (|9.10H9.17|) agree with the results produced by the integrals in Eqs. (j4.8M.13j) as 
listed in Eas. (l9.1H9.7l) . 
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